3.1.14 \(\int \frac {(A+B x) (a+b x+c x^2)^2}{d+e x+f x^2} \, dx\) [14]

3.1.14.1 Optimal result
3.1.14.2 Mathematica [A] (verified)
3.1.14.3 Rubi [A] (verified)
3.1.14.4 Maple [A] (verified)
3.1.14.5 Fricas [A] (verification not implemented)
3.1.14.6 Sympy [B] (verification not implemented)
3.1.14.7 Maxima [F(-2)]
3.1.14.8 Giac [A] (verification not implemented)
3.1.14.9 Mupad [B] (verification not implemented)

3.1.14.1 Optimal result

Integrand size = 30, antiderivative size = 542 \[ \int \frac {(A+B x) \left (a+b x+c x^2\right )^2}{d+e x+f x^2} \, dx=\frac {\left (B (c e-b f) \left (f (b e-2 a f)-c \left (e^2-2 d f\right )\right )+A f \left (b^2 f^2-2 c f (b e-a f)+c^2 \left (e^2-d f\right )\right )\right ) x}{f^4}-\frac {\left (A c f (c e-2 b f)-B \left (b^2 f^2-2 c f (b e-a f)+c^2 \left (e^2-d f\right )\right )\right ) x^2}{2 f^3}-\frac {c (B c e-2 b B f-A c f) x^3}{3 f^2}+\frac {B c^2 x^4}{4 f}-\frac {\left (A f \left (c^2 \left (e^4-4 d e^2 f+2 d^2 f^2\right )-f^2 \left (2 a b e f-2 a^2 f^2-b^2 \left (e^2-2 d f\right )\right )+2 c f \left (a f \left (e^2-2 d f\right )-b \left (e^3-3 d e f\right )\right )\right )-B \left (c^2 \left (e^5-5 d e^3 f+5 d^2 e f^2\right )+f^2 \left (a^2 e f^2-2 a b f \left (e^2-2 d f\right )+b^2 \left (e^3-3 d e f\right )\right )+2 c f \left (a e f \left (e^2-3 d f\right )-b \left (e^4-4 d e^2 f+2 d^2 f^2\right )\right )\right )\right ) \text {arctanh}\left (\frac {e+2 f x}{\sqrt {e^2-4 d f}}\right )}{f^5 \sqrt {e^2-4 d f}}+\frac {\left (A f (c e-b f) \left (f (b e-2 a f)-c \left (e^2-2 d f\right )\right )+B \left (c^2 \left (e^4-3 d e^2 f+d^2 f^2\right )-f^2 \left (2 a b e f-a^2 f^2-b^2 \left (e^2-d f\right )\right )+2 c f \left (a f \left (e^2-d f\right )-b \left (e^3-2 d e f\right )\right )\right )\right ) \log \left (d+e x+f x^2\right )}{2 f^5} \]

output
(B*(-b*f+c*e)*(f*(-2*a*f+b*e)-c*(-2*d*f+e^2))+A*f*(b^2*f^2-2*c*f*(-a*f+b*e 
)+c^2*(-d*f+e^2)))*x/f^4-1/2*(A*c*f*(-2*b*f+c*e)-B*(b^2*f^2-2*c*f*(-a*f+b* 
e)+c^2*(-d*f+e^2)))*x^2/f^3-1/3*c*(-A*c*f-2*B*b*f+B*c*e)*x^3/f^2+1/4*B*c^2 
*x^4/f+1/2*(A*f*(-b*f+c*e)*(f*(-2*a*f+b*e)-c*(-2*d*f+e^2))+B*(c^2*(d^2*f^2 
-3*d*e^2*f+e^4)-f^2*(2*a*b*e*f-a^2*f^2-b^2*(-d*f+e^2))+2*c*f*(a*f*(-d*f+e^ 
2)-b*(-2*d*e*f+e^3))))*ln(f*x^2+e*x+d)/f^5-(A*f*(c^2*(2*d^2*f^2-4*d*e^2*f+ 
e^4)-f^2*(2*a*b*e*f-2*a^2*f^2-b^2*(-2*d*f+e^2))+2*c*f*(a*f*(-2*d*f+e^2)-b* 
(-3*d*e*f+e^3)))-B*(c^2*(5*d^2*e*f^2-5*d*e^3*f+e^5)+f^2*(a^2*e*f^2-2*a*b*f 
*(-2*d*f+e^2)+b^2*(-3*d*e*f+e^3))+2*c*f*(a*e*f*(-3*d*f+e^2)-b*(2*d^2*f^2-4 
*d*e^2*f+e^4))))*arctanh((2*f*x+e)/(-4*d*f+e^2)^(1/2))/f^5/(-4*d*f+e^2)^(1 
/2)
 
3.1.14.2 Mathematica [A] (verified)

Time = 0.45 (sec) , antiderivative size = 535, normalized size of antiderivative = 0.99 \[ \int \frac {(A+B x) \left (a+b x+c x^2\right )^2}{d+e x+f x^2} \, dx=\frac {12 f \left (-B (c e-b f) \left (f (-b e+2 a f)+c \left (e^2-2 d f\right )\right )+A f \left (b^2 f^2+2 c f (-b e+a f)+c^2 \left (e^2-d f\right )\right )\right ) x+6 f^2 \left (A c f (-c e+2 b f)+B \left (b^2 f^2+2 c f (-b e+a f)+c^2 \left (e^2-d f\right )\right )\right ) x^2+4 c f^3 (-B c e+2 b B f+A c f) x^3+3 B c^2 f^4 x^4-\frac {12 \left (-A f \left (c^2 \left (e^4-4 d e^2 f+2 d^2 f^2\right )+f^2 \left (-2 a b e f+2 a^2 f^2+b^2 \left (e^2-2 d f\right )\right )+2 c f \left (a f \left (e^2-2 d f\right )-b \left (e^3-3 d e f\right )\right )\right )+B \left (c^2 \left (e^5-5 d e^3 f+5 d^2 e f^2\right )+f^2 \left (a^2 e f^2+2 a b f \left (-e^2+2 d f\right )+b^2 \left (e^3-3 d e f\right )\right )-2 c f \left (-a e f \left (e^2-3 d f\right )+b \left (e^4-4 d e^2 f+2 d^2 f^2\right )\right )\right )\right ) \arctan \left (\frac {e+2 f x}{\sqrt {-e^2+4 d f}}\right )}{\sqrt {-e^2+4 d f}}+6 \left (A f (-c e+b f) \left (f (-b e+2 a f)+c \left (e^2-2 d f\right )\right )+B \left (c^2 \left (e^4-3 d e^2 f+d^2 f^2\right )+f^2 \left (-2 a b e f+a^2 f^2+b^2 \left (e^2-d f\right )\right )-2 c f \left (a f \left (-e^2+d f\right )+b \left (e^3-2 d e f\right )\right )\right )\right ) \log (d+x (e+f x))}{12 f^5} \]

input
Integrate[((A + B*x)*(a + b*x + c*x^2)^2)/(d + e*x + f*x^2),x]
 
output
(12*f*(-(B*(c*e - b*f)*(f*(-(b*e) + 2*a*f) + c*(e^2 - 2*d*f))) + A*f*(b^2* 
f^2 + 2*c*f*(-(b*e) + a*f) + c^2*(e^2 - d*f)))*x + 6*f^2*(A*c*f*(-(c*e) + 
2*b*f) + B*(b^2*f^2 + 2*c*f*(-(b*e) + a*f) + c^2*(e^2 - d*f)))*x^2 + 4*c*f 
^3*(-(B*c*e) + 2*b*B*f + A*c*f)*x^3 + 3*B*c^2*f^4*x^4 - (12*(-(A*f*(c^2*(e 
^4 - 4*d*e^2*f + 2*d^2*f^2) + f^2*(-2*a*b*e*f + 2*a^2*f^2 + b^2*(e^2 - 2*d 
*f)) + 2*c*f*(a*f*(e^2 - 2*d*f) - b*(e^3 - 3*d*e*f)))) + B*(c^2*(e^5 - 5*d 
*e^3*f + 5*d^2*e*f^2) + f^2*(a^2*e*f^2 + 2*a*b*f*(-e^2 + 2*d*f) + b^2*(e^3 
 - 3*d*e*f)) - 2*c*f*(-(a*e*f*(e^2 - 3*d*f)) + b*(e^4 - 4*d*e^2*f + 2*d^2* 
f^2))))*ArcTan[(e + 2*f*x)/Sqrt[-e^2 + 4*d*f]])/Sqrt[-e^2 + 4*d*f] + 6*(A* 
f*(-(c*e) + b*f)*(f*(-(b*e) + 2*a*f) + c*(e^2 - 2*d*f)) + B*(c^2*(e^4 - 3* 
d*e^2*f + d^2*f^2) + f^2*(-2*a*b*e*f + a^2*f^2 + b^2*(e^2 - d*f)) - 2*c*f* 
(a*f*(-e^2 + d*f) + b*(e^3 - 2*d*e*f))))*Log[d + x*(e + f*x)])/(12*f^5)
 
3.1.14.3 Rubi [A] (verified)

Time = 1.11 (sec) , antiderivative size = 542, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.067, Rules used = {1344, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {(A+B x) \left (a+b x+c x^2\right )^2}{d+e x+f x^2} \, dx\)

\(\Big \downarrow \) 1344

\(\displaystyle \int \left (\frac {x \left (B \left (-f^2 \left (-a^2 f^2+2 a b e f-\left (b^2 \left (e^2-d f\right )\right )\right )+2 c f \left (a f \left (e^2-d f\right )-b \left (e^3-2 d e f\right )\right )+c^2 \left (d^2 f^2-3 d e^2 f+e^4\right )\right )+A f (c e-b f) \left (f (b e-2 a f)-c \left (e^2-2 d f\right )\right )\right )+A f \left (-f^2 \left (b^2 d-a^2 f\right )+2 c d f (b e-a f)+c^2 (-d) \left (e^2-d f\right )\right )-B d (c e-b f) \left (f (b e-2 a f)-c \left (e^2-2 d f\right )\right )}{f^4 \left (d+e x+f x^2\right )}+\frac {A f \left (-2 c f (b e-a f)+b^2 f^2+c^2 \left (e^2-d f\right )\right )+B (c e-b f) \left (f (b e-2 a f)-c \left (e^2-2 d f\right )\right )}{f^4}-\frac {x \left (A c f (c e-2 b f)-B \left (-2 c f (b e-a f)+b^2 f^2+c^2 \left (e^2-d f\right )\right )\right )}{f^3}-\frac {c x^2 (-A c f-2 b B f+B c e)}{f^2}+\frac {B c^2 x^3}{f}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle -\frac {\text {arctanh}\left (\frac {e+2 f x}{\sqrt {e^2-4 d f}}\right ) \left (A f \left (-f^2 \left (-2 a^2 f^2+2 a b e f-\left (b^2 \left (e^2-2 d f\right )\right )\right )+2 c f \left (a f \left (e^2-2 d f\right )-b \left (e^3-3 d e f\right )\right )+c^2 \left (2 d^2 f^2-4 d e^2 f+e^4\right )\right )-B \left (f^2 \left (a^2 e f^2-2 a b f \left (e^2-2 d f\right )+b^2 \left (e^3-3 d e f\right )\right )+2 c f \left (a e f \left (e^2-3 d f\right )-b \left (2 d^2 f^2-4 d e^2 f+e^4\right )\right )+c^2 \left (5 d^2 e f^2-5 d e^3 f+e^5\right )\right )\right )}{f^5 \sqrt {e^2-4 d f}}+\frac {\log \left (d+e x+f x^2\right ) \left (B \left (-f^2 \left (-a^2 f^2+2 a b e f-\left (b^2 \left (e^2-d f\right )\right )\right )+2 c f \left (a f \left (e^2-d f\right )-b \left (e^3-2 d e f\right )\right )+c^2 \left (d^2 f^2-3 d e^2 f+e^4\right )\right )+A f (c e-b f) \left (f (b e-2 a f)-c \left (e^2-2 d f\right )\right )\right )}{2 f^5}+\frac {x \left (A f \left (-2 c f (b e-a f)+b^2 f^2+c^2 \left (e^2-d f\right )\right )+B (c e-b f) \left (f (b e-2 a f)-c \left (e^2-2 d f\right )\right )\right )}{f^4}-\frac {x^2 \left (A c f (c e-2 b f)-B \left (-2 c f (b e-a f)+b^2 f^2+c^2 \left (e^2-d f\right )\right )\right )}{2 f^3}-\frac {c x^3 (-A c f-2 b B f+B c e)}{3 f^2}+\frac {B c^2 x^4}{4 f}\)

input
Int[((A + B*x)*(a + b*x + c*x^2)^2)/(d + e*x + f*x^2),x]
 
output
((B*(c*e - b*f)*(f*(b*e - 2*a*f) - c*(e^2 - 2*d*f)) + A*f*(b^2*f^2 - 2*c*f 
*(b*e - a*f) + c^2*(e^2 - d*f)))*x)/f^4 - ((A*c*f*(c*e - 2*b*f) - B*(b^2*f 
^2 - 2*c*f*(b*e - a*f) + c^2*(e^2 - d*f)))*x^2)/(2*f^3) - (c*(B*c*e - 2*b* 
B*f - A*c*f)*x^3)/(3*f^2) + (B*c^2*x^4)/(4*f) - ((A*f*(c^2*(e^4 - 4*d*e^2* 
f + 2*d^2*f^2) - f^2*(2*a*b*e*f - 2*a^2*f^2 - b^2*(e^2 - 2*d*f)) + 2*c*f*( 
a*f*(e^2 - 2*d*f) - b*(e^3 - 3*d*e*f))) - B*(c^2*(e^5 - 5*d*e^3*f + 5*d^2* 
e*f^2) + f^2*(a^2*e*f^2 - 2*a*b*f*(e^2 - 2*d*f) + b^2*(e^3 - 3*d*e*f)) + 2 
*c*f*(a*e*f*(e^2 - 3*d*f) - b*(e^4 - 4*d*e^2*f + 2*d^2*f^2))))*ArcTanh[(e 
+ 2*f*x)/Sqrt[e^2 - 4*d*f]])/(f^5*Sqrt[e^2 - 4*d*f]) + ((A*f*(c*e - b*f)*( 
f*(b*e - 2*a*f) - c*(e^2 - 2*d*f)) + B*(c^2*(e^4 - 3*d*e^2*f + d^2*f^2) - 
f^2*(2*a*b*e*f - a^2*f^2 - b^2*(e^2 - d*f)) + 2*c*f*(a*f*(e^2 - d*f) - b*( 
e^3 - 2*d*e*f))))*Log[d + e*x + f*x^2])/(2*f^5)
 

3.1.14.3.1 Defintions of rubi rules used

rule 1344
Int[((g_.) + (h_.)*(x_))*((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_)*((d_) + (e 
_.)*(x_) + (f_.)*(x_)^2)^(q_), x_Symbol] :> Int[ExpandIntegrand[(a + b*x + 
c*x^2)^p*(d + e*x + f*x^2)^q*(g + h*x), x], x] /; FreeQ[{a, b, c, d, e, f, 
g, h}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[e^2 - 4*d*f, 0] && IGtQ[p, 0] && In 
tegerQ[q]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.1.14.4 Maple [A] (verified)

Time = 0.92 (sec) , antiderivative size = 800, normalized size of antiderivative = 1.48

method result size
default \(\frac {\frac {1}{4} B \,c^{2} x^{4} f^{3}+\frac {1}{3} A \,c^{2} f^{3} x^{3}+\frac {2}{3} B b c \,f^{3} x^{3}-\frac {1}{3} B \,c^{2} e \,f^{2} x^{3}+A b c \,f^{3} x^{2}-\frac {1}{2} A \,c^{2} e \,f^{2} x^{2}+B a c \,f^{3} x^{2}+\frac {1}{2} B \,b^{2} f^{3} x^{2}-B b c e \,f^{2} x^{2}-\frac {1}{2} B \,c^{2} d \,f^{2} x^{2}+\frac {1}{2} B \,c^{2} e^{2} f \,x^{2}+2 A a c \,f^{3} x +A \,b^{2} f^{3} x -2 A b c e \,f^{2} x -A \,c^{2} d \,f^{2} x +A \,c^{2} e^{2} f x +2 B a b \,f^{3} x -2 B a c e \,f^{2} x -B \,b^{2} e \,f^{2} x -2 B b c d \,f^{2} x +2 B b c \,e^{2} f x +2 B \,c^{2} d e f x -B \,c^{2} e^{3} x}{f^{4}}+\frac {\frac {\left (2 A a b \,f^{4}-2 A a c e \,f^{3}-A \,b^{2} e \,f^{3}-2 A b c d \,f^{3}+2 A b c \,e^{2} f^{2}+2 A \,c^{2} d e \,f^{2}-A \,c^{2} e^{3} f +B \,a^{2} f^{4}-2 B a b e \,f^{3}-2 B a c d \,f^{3}+2 B a c \,e^{2} f^{2}-B \,b^{2} d \,f^{3}+B \,b^{2} e^{2} f^{2}+4 B b c d e \,f^{2}-2 B b c \,e^{3} f +B \,c^{2} d^{2} f^{2}-3 B \,c^{2} d \,e^{2} f +B \,c^{2} e^{4}\right ) \ln \left (f \,x^{2}+e x +d \right )}{2 f}+\frac {2 \left (A \,a^{2} f^{4}-2 A a c d \,f^{3}-A \,b^{2} d \,f^{3}+2 A b c d e \,f^{2}+A \,c^{2} d^{2} f^{2}-A \,c^{2} d \,e^{2} f -2 B a b d \,f^{3}+2 B a c d e \,f^{2}+B \,b^{2} d e \,f^{2}+2 B b c \,d^{2} f^{2}-2 B b c d \,e^{2} f -2 B \,c^{2} d^{2} e f +B \,c^{2} d \,e^{3}-\frac {\left (2 A a b \,f^{4}-2 A a c e \,f^{3}-A \,b^{2} e \,f^{3}-2 A b c d \,f^{3}+2 A b c \,e^{2} f^{2}+2 A \,c^{2} d e \,f^{2}-A \,c^{2} e^{3} f +B \,a^{2} f^{4}-2 B a b e \,f^{3}-2 B a c d \,f^{3}+2 B a c \,e^{2} f^{2}-B \,b^{2} d \,f^{3}+B \,b^{2} e^{2} f^{2}+4 B b c d e \,f^{2}-2 B b c \,e^{3} f +B \,c^{2} d^{2} f^{2}-3 B \,c^{2} d \,e^{2} f +B \,c^{2} e^{4}\right ) e}{2 f}\right ) \arctan \left (\frac {2 f x +e}{\sqrt {4 d f -e^{2}}}\right )}{\sqrt {4 d f -e^{2}}}}{f^{4}}\) \(800\)
risch \(\text {Expression too large to display}\) \(60299\)

input
int((B*x+A)*(c*x^2+b*x+a)^2/(f*x^2+e*x+d),x,method=_RETURNVERBOSE)
 
output
1/f^4*(1/4*B*c^2*x^4*f^3+1/3*A*c^2*f^3*x^3+2/3*B*b*c*f^3*x^3-1/3*B*c^2*e*f 
^2*x^3+A*b*c*f^3*x^2-1/2*A*c^2*e*f^2*x^2+B*a*c*f^3*x^2+1/2*B*b^2*f^3*x^2-B 
*b*c*e*f^2*x^2-1/2*B*c^2*d*f^2*x^2+1/2*B*c^2*e^2*f*x^2+2*A*a*c*f^3*x+A*b^2 
*f^3*x-2*A*b*c*e*f^2*x-A*c^2*d*f^2*x+A*c^2*e^2*f*x+2*B*a*b*f^3*x-2*B*a*c*e 
*f^2*x-B*b^2*e*f^2*x-2*B*b*c*d*f^2*x+2*B*b*c*e^2*f*x+2*B*c^2*d*e*f*x-B*c^2 
*e^3*x)+1/f^4*(1/2*(2*A*a*b*f^4-2*A*a*c*e*f^3-A*b^2*e*f^3-2*A*b*c*d*f^3+2* 
A*b*c*e^2*f^2+2*A*c^2*d*e*f^2-A*c^2*e^3*f+B*a^2*f^4-2*B*a*b*e*f^3-2*B*a*c* 
d*f^3+2*B*a*c*e^2*f^2-B*b^2*d*f^3+B*b^2*e^2*f^2+4*B*b*c*d*e*f^2-2*B*b*c*e^ 
3*f+B*c^2*d^2*f^2-3*B*c^2*d*e^2*f+B*c^2*e^4)/f*ln(f*x^2+e*x+d)+2*(A*a^2*f^ 
4-2*A*a*c*d*f^3-A*b^2*d*f^3+2*A*b*c*d*e*f^2+A*c^2*d^2*f^2-A*c^2*d*e^2*f-2* 
B*a*b*d*f^3+2*B*a*c*d*e*f^2+B*b^2*d*e*f^2+2*B*b*c*d^2*f^2-2*B*b*c*d*e^2*f- 
2*B*c^2*d^2*e*f+B*c^2*d*e^3-1/2*(2*A*a*b*f^4-2*A*a*c*e*f^3-A*b^2*e*f^3-2*A 
*b*c*d*f^3+2*A*b*c*e^2*f^2+2*A*c^2*d*e*f^2-A*c^2*e^3*f+B*a^2*f^4-2*B*a*b*e 
*f^3-2*B*a*c*d*f^3+2*B*a*c*e^2*f^2-B*b^2*d*f^3+B*b^2*e^2*f^2+4*B*b*c*d*e*f 
^2-2*B*b*c*e^3*f+B*c^2*d^2*f^2-3*B*c^2*d*e^2*f+B*c^2*e^4)*e/f)/(4*d*f-e^2) 
^(1/2)*arctan((2*f*x+e)/(4*d*f-e^2)^(1/2)))
 
3.1.14.5 Fricas [A] (verification not implemented)

Time = 0.42 (sec) , antiderivative size = 1837, normalized size of antiderivative = 3.39 \[ \int \frac {(A+B x) \left (a+b x+c x^2\right )^2}{d+e x+f x^2} \, dx=\text {Too large to display} \]

input
integrate((B*x+A)*(c*x^2+b*x+a)^2/(f*x^2+e*x+d),x, algorithm="fricas")
 
output
[1/12*(3*(B*c^2*e^2*f^4 - 4*B*c^2*d*f^5)*x^4 - 4*(B*c^2*e^3*f^3 + 4*(2*B*b 
*c + A*c^2)*d*f^5 - (4*B*c^2*d*e + (2*B*b*c + A*c^2)*e^2)*f^4)*x^3 + 6*(B* 
c^2*e^4*f^2 - 4*(B*b^2 + 2*(B*a + A*b)*c)*d*f^5 + (4*B*c^2*d^2 + 4*(2*B*b* 
c + A*c^2)*d*e + (B*b^2 + 2*(B*a + A*b)*c)*e^2)*f^4 - (5*B*c^2*d*e^2 + (2* 
B*b*c + A*c^2)*e^3)*f^3)*x^2 - 6*(B*c^2*e^5 - 2*A*a^2*f^5 + (2*(2*B*a*b + 
A*b^2 + 2*A*a*c)*d + (B*a^2 + 2*A*a*b)*e)*f^4 - (2*(2*B*b*c + A*c^2)*d^2 + 
 3*(B*b^2 + 2*(B*a + A*b)*c)*d*e + (2*B*a*b + A*b^2 + 2*A*a*c)*e^2)*f^3 + 
(5*B*c^2*d^2*e + 4*(2*B*b*c + A*c^2)*d*e^2 + (B*b^2 + 2*(B*a + A*b)*c)*e^3 
)*f^2 - (5*B*c^2*d*e^3 + (2*B*b*c + A*c^2)*e^4)*f)*sqrt(e^2 - 4*d*f)*log(( 
2*f^2*x^2 + 2*e*f*x + e^2 - 2*d*f - sqrt(e^2 - 4*d*f)*(2*f*x + e))/(f*x^2 
+ e*x + d)) - 12*(B*c^2*e^5*f + 4*(2*B*a*b + A*b^2 + 2*A*a*c)*d*f^5 - (4*( 
2*B*b*c + A*c^2)*d^2 + 4*(B*b^2 + 2*(B*a + A*b)*c)*d*e + (2*B*a*b + A*b^2 
+ 2*A*a*c)*e^2)*f^4 + (8*B*c^2*d^2*e + 5*(2*B*b*c + A*c^2)*d*e^2 + (B*b^2 
+ 2*(B*a + A*b)*c)*e^3)*f^3 - (6*B*c^2*d*e^3 + (2*B*b*c + A*c^2)*e^4)*f^2) 
*x + 6*(B*c^2*e^6 - 4*(B*a^2 + 2*A*a*b)*d*f^5 + (4*(B*b^2 + 2*(B*a + A*b)* 
c)*d^2 + 4*(2*B*a*b + A*b^2 + 2*A*a*c)*d*e + (B*a^2 + 2*A*a*b)*e^2)*f^4 - 
(4*B*c^2*d^3 + 8*(2*B*b*c + A*c^2)*d^2*e + 5*(B*b^2 + 2*(B*a + A*b)*c)*d*e 
^2 + (2*B*a*b + A*b^2 + 2*A*a*c)*e^3)*f^3 + (13*B*c^2*d^2*e^2 + 6*(2*B*b*c 
 + A*c^2)*d*e^3 + (B*b^2 + 2*(B*a + A*b)*c)*e^4)*f^2 - (7*B*c^2*d*e^4 + (2 
*B*b*c + A*c^2)*e^5)*f)*log(f*x^2 + e*x + d))/(e^2*f^5 - 4*d*f^6), 1/12...
 
3.1.14.6 Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 4663 vs. \(2 (520) = 1040\).

Time = 82.39 (sec) , antiderivative size = 4663, normalized size of antiderivative = 8.60 \[ \int \frac {(A+B x) \left (a+b x+c x^2\right )^2}{d+e x+f x^2} \, dx=\text {Too large to display} \]

input
integrate((B*x+A)*(c*x**2+b*x+a)**2/(f*x**2+e*x+d),x)
 
output
B*c**2*x**4/(4*f) + x**3*(A*c**2/(3*f) + 2*B*b*c/(3*f) - B*c**2*e/(3*f**2) 
) + x**2*(A*b*c/f - A*c**2*e/(2*f**2) + B*a*c/f + B*b**2/(2*f) - B*b*c*e/f 
**2 - B*c**2*d/(2*f**2) + B*c**2*e**2/(2*f**3)) + x*(2*A*a*c/f + A*b**2/f 
- 2*A*b*c*e/f**2 - A*c**2*d/f**2 + A*c**2*e**2/f**3 + 2*B*a*b/f - 2*B*a*c* 
e/f**2 - B*b**2*e/f**2 - 2*B*b*c*d/f**2 + 2*B*b*c*e**2/f**3 + 2*B*c**2*d*e 
/f**3 - B*c**2*e**3/f**4) + (-sqrt(-4*d*f + e**2)*(-2*A*a**2*f**5 + 2*A*a* 
b*e*f**4 + 4*A*a*c*d*f**4 - 2*A*a*c*e**2*f**3 + 2*A*b**2*d*f**4 - A*b**2*e 
**2*f**3 - 6*A*b*c*d*e*f**3 + 2*A*b*c*e**3*f**2 - 2*A*c**2*d**2*f**3 + 4*A 
*c**2*d*e**2*f**2 - A*c**2*e**4*f + B*a**2*e*f**4 + 4*B*a*b*d*f**4 - 2*B*a 
*b*e**2*f**3 - 6*B*a*c*d*e*f**3 + 2*B*a*c*e**3*f**2 - 3*B*b**2*d*e*f**3 + 
B*b**2*e**3*f**2 - 4*B*b*c*d**2*f**3 + 8*B*b*c*d*e**2*f**2 - 2*B*b*c*e**4* 
f + 5*B*c**2*d**2*e*f**2 - 5*B*c**2*d*e**3*f + B*c**2*e**5)/(2*f**5*(4*d*f 
 - e**2)) + (2*A*a*b*f**4 - 2*A*a*c*e*f**3 - A*b**2*e*f**3 - 2*A*b*c*d*f** 
3 + 2*A*b*c*e**2*f**2 + 2*A*c**2*d*e*f**2 - A*c**2*e**3*f + B*a**2*f**4 - 
2*B*a*b*e*f**3 - 2*B*a*c*d*f**3 + 2*B*a*c*e**2*f**2 - B*b**2*d*f**3 + B*b* 
*2*e**2*f**2 + 4*B*b*c*d*e*f**2 - 2*B*b*c*e**3*f + B*c**2*d**2*f**2 - 3*B* 
c**2*d*e**2*f + B*c**2*e**4)/(2*f**5))*log(x + (-A*a**2*e*f**4 + 4*A*a*b*d 
*f**4 - 2*A*a*c*d*e*f**3 - A*b**2*d*e*f**3 - 4*A*b*c*d**2*f**3 + 2*A*b*c*d 
*e**2*f**2 + 3*A*c**2*d**2*e*f**2 - A*c**2*d*e**3*f + 2*B*a**2*d*f**4 - 2* 
B*a*b*d*e*f**3 - 4*B*a*c*d**2*f**3 + 2*B*a*c*d*e**2*f**2 - 2*B*b**2*d**...
 
3.1.14.7 Maxima [F(-2)]

Exception generated. \[ \int \frac {(A+B x) \left (a+b x+c x^2\right )^2}{d+e x+f x^2} \, dx=\text {Exception raised: ValueError} \]

input
integrate((B*x+A)*(c*x^2+b*x+a)^2/(f*x^2+e*x+d),x, algorithm="maxima")
 
output
Exception raised: ValueError >> Computation failed since Maxima requested 
additional constraints; using the 'assume' command before evaluation *may* 
 help (example of legal syntax is 'assume(4*d*f-e^2>0)', see `assume?` for 
 more deta
 
3.1.14.8 Giac [A] (verification not implemented)

Time = 0.28 (sec) , antiderivative size = 743, normalized size of antiderivative = 1.37 \[ \int \frac {(A+B x) \left (a+b x+c x^2\right )^2}{d+e x+f x^2} \, dx=\frac {3 \, B c^{2} f^{3} x^{4} - 4 \, B c^{2} e f^{2} x^{3} + 8 \, B b c f^{3} x^{3} + 4 \, A c^{2} f^{3} x^{3} + 6 \, B c^{2} e^{2} f x^{2} - 6 \, B c^{2} d f^{2} x^{2} - 12 \, B b c e f^{2} x^{2} - 6 \, A c^{2} e f^{2} x^{2} + 6 \, B b^{2} f^{3} x^{2} + 12 \, B a c f^{3} x^{2} + 12 \, A b c f^{3} x^{2} - 12 \, B c^{2} e^{3} x + 24 \, B c^{2} d e f x + 24 \, B b c e^{2} f x + 12 \, A c^{2} e^{2} f x - 24 \, B b c d f^{2} x - 12 \, A c^{2} d f^{2} x - 12 \, B b^{2} e f^{2} x - 24 \, B a c e f^{2} x - 24 \, A b c e f^{2} x + 24 \, B a b f^{3} x + 12 \, A b^{2} f^{3} x + 24 \, A a c f^{3} x}{12 \, f^{4}} + \frac {{\left (B c^{2} e^{4} - 3 \, B c^{2} d e^{2} f - 2 \, B b c e^{3} f - A c^{2} e^{3} f + B c^{2} d^{2} f^{2} + 4 \, B b c d e f^{2} + 2 \, A c^{2} d e f^{2} + B b^{2} e^{2} f^{2} + 2 \, B a c e^{2} f^{2} + 2 \, A b c e^{2} f^{2} - B b^{2} d f^{3} - 2 \, B a c d f^{3} - 2 \, A b c d f^{3} - 2 \, B a b e f^{3} - A b^{2} e f^{3} - 2 \, A a c e f^{3} + B a^{2} f^{4} + 2 \, A a b f^{4}\right )} \log \left (f x^{2} + e x + d\right )}{2 \, f^{5}} - \frac {{\left (B c^{2} e^{5} - 5 \, B c^{2} d e^{3} f - 2 \, B b c e^{4} f - A c^{2} e^{4} f + 5 \, B c^{2} d^{2} e f^{2} + 8 \, B b c d e^{2} f^{2} + 4 \, A c^{2} d e^{2} f^{2} + B b^{2} e^{3} f^{2} + 2 \, B a c e^{3} f^{2} + 2 \, A b c e^{3} f^{2} - 4 \, B b c d^{2} f^{3} - 2 \, A c^{2} d^{2} f^{3} - 3 \, B b^{2} d e f^{3} - 6 \, B a c d e f^{3} - 6 \, A b c d e f^{3} - 2 \, B a b e^{2} f^{3} - A b^{2} e^{2} f^{3} - 2 \, A a c e^{2} f^{3} + 4 \, B a b d f^{4} + 2 \, A b^{2} d f^{4} + 4 \, A a c d f^{4} + B a^{2} e f^{4} + 2 \, A a b e f^{4} - 2 \, A a^{2} f^{5}\right )} \arctan \left (\frac {2 \, f x + e}{\sqrt {-e^{2} + 4 \, d f}}\right )}{\sqrt {-e^{2} + 4 \, d f} f^{5}} \]

input
integrate((B*x+A)*(c*x^2+b*x+a)^2/(f*x^2+e*x+d),x, algorithm="giac")
 
output
1/12*(3*B*c^2*f^3*x^4 - 4*B*c^2*e*f^2*x^3 + 8*B*b*c*f^3*x^3 + 4*A*c^2*f^3* 
x^3 + 6*B*c^2*e^2*f*x^2 - 6*B*c^2*d*f^2*x^2 - 12*B*b*c*e*f^2*x^2 - 6*A*c^2 
*e*f^2*x^2 + 6*B*b^2*f^3*x^2 + 12*B*a*c*f^3*x^2 + 12*A*b*c*f^3*x^2 - 12*B* 
c^2*e^3*x + 24*B*c^2*d*e*f*x + 24*B*b*c*e^2*f*x + 12*A*c^2*e^2*f*x - 24*B* 
b*c*d*f^2*x - 12*A*c^2*d*f^2*x - 12*B*b^2*e*f^2*x - 24*B*a*c*e*f^2*x - 24* 
A*b*c*e*f^2*x + 24*B*a*b*f^3*x + 12*A*b^2*f^3*x + 24*A*a*c*f^3*x)/f^4 + 1/ 
2*(B*c^2*e^4 - 3*B*c^2*d*e^2*f - 2*B*b*c*e^3*f - A*c^2*e^3*f + B*c^2*d^2*f 
^2 + 4*B*b*c*d*e*f^2 + 2*A*c^2*d*e*f^2 + B*b^2*e^2*f^2 + 2*B*a*c*e^2*f^2 + 
 2*A*b*c*e^2*f^2 - B*b^2*d*f^3 - 2*B*a*c*d*f^3 - 2*A*b*c*d*f^3 - 2*B*a*b*e 
*f^3 - A*b^2*e*f^3 - 2*A*a*c*e*f^3 + B*a^2*f^4 + 2*A*a*b*f^4)*log(f*x^2 + 
e*x + d)/f^5 - (B*c^2*e^5 - 5*B*c^2*d*e^3*f - 2*B*b*c*e^4*f - A*c^2*e^4*f 
+ 5*B*c^2*d^2*e*f^2 + 8*B*b*c*d*e^2*f^2 + 4*A*c^2*d*e^2*f^2 + B*b^2*e^3*f^ 
2 + 2*B*a*c*e^3*f^2 + 2*A*b*c*e^3*f^2 - 4*B*b*c*d^2*f^3 - 2*A*c^2*d^2*f^3 
- 3*B*b^2*d*e*f^3 - 6*B*a*c*d*e*f^3 - 6*A*b*c*d*e*f^3 - 2*B*a*b*e^2*f^3 - 
A*b^2*e^2*f^3 - 2*A*a*c*e^2*f^3 + 4*B*a*b*d*f^4 + 2*A*b^2*d*f^4 + 4*A*a*c* 
d*f^4 + B*a^2*e*f^4 + 2*A*a*b*e*f^4 - 2*A*a^2*f^5)*arctan((2*f*x + e)/sqrt 
(-e^2 + 4*d*f))/(sqrt(-e^2 + 4*d*f)*f^5)
 
3.1.14.9 Mupad [B] (verification not implemented)

Time = 15.37 (sec) , antiderivative size = 893, normalized size of antiderivative = 1.65 \[ \int \frac {(A+B x) \left (a+b x+c x^2\right )^2}{d+e x+f x^2} \, dx=x^3\,\left (\frac {A\,c^2+2\,B\,b\,c}{3\,f}-\frac {B\,c^2\,e}{3\,f^2}\right )+x\,\left (\frac {A\,b^2+2\,B\,a\,b+2\,A\,a\,c}{f}-\frac {d\,\left (\frac {A\,c^2+2\,B\,b\,c}{f}-\frac {B\,c^2\,e}{f^2}\right )}{f}+\frac {e\,\left (\frac {e\,\left (\frac {A\,c^2+2\,B\,b\,c}{f}-\frac {B\,c^2\,e}{f^2}\right )}{f}-\frac {B\,b^2+2\,A\,c\,b+2\,B\,a\,c}{f}+\frac {B\,c^2\,d}{f^2}\right )}{f}\right )-x^2\,\left (\frac {e\,\left (\frac {A\,c^2+2\,B\,b\,c}{f}-\frac {B\,c^2\,e}{f^2}\right )}{2\,f}-\frac {B\,b^2+2\,A\,c\,b+2\,B\,a\,c}{2\,f}+\frac {B\,c^2\,d}{2\,f^2}\right )-\frac {\ln \left (f\,x^2+e\,x+d\right )\,\left (-4\,B\,a^2\,d\,f^5+B\,a^2\,e^2\,f^4+8\,B\,a\,b\,d\,e\,f^4-8\,A\,a\,b\,d\,f^5-2\,B\,a\,b\,e^3\,f^3+2\,A\,a\,b\,e^2\,f^4+8\,B\,a\,c\,d^2\,f^4-10\,B\,a\,c\,d\,e^2\,f^3+8\,A\,a\,c\,d\,e\,f^4+2\,B\,a\,c\,e^4\,f^2-2\,A\,a\,c\,e^3\,f^3+4\,B\,b^2\,d^2\,f^4-5\,B\,b^2\,d\,e^2\,f^3+4\,A\,b^2\,d\,e\,f^4+B\,b^2\,e^4\,f^2-A\,b^2\,e^3\,f^3-16\,B\,b\,c\,d^2\,e\,f^3+8\,A\,b\,c\,d^2\,f^4+12\,B\,b\,c\,d\,e^3\,f^2-10\,A\,b\,c\,d\,e^2\,f^3-2\,B\,b\,c\,e^5\,f+2\,A\,b\,c\,e^4\,f^2-4\,B\,c^2\,d^3\,f^3+13\,B\,c^2\,d^2\,e^2\,f^2-8\,A\,c^2\,d^2\,e\,f^3-7\,B\,c^2\,d\,e^4\,f+6\,A\,c^2\,d\,e^3\,f^2+B\,c^2\,e^6-A\,c^2\,e^5\,f\right )}{2\,\left (4\,d\,f^6-e^2\,f^5\right )}+\frac {B\,c^2\,x^4}{4\,f}+\frac {\mathrm {atan}\left (\frac {e}{\sqrt {4\,d\,f-e^2}}+\frac {2\,f\,x}{\sqrt {4\,d\,f-e^2}}\right )\,\left (-B\,a^2\,e\,f^4+2\,A\,a^2\,f^5-4\,B\,a\,b\,d\,f^4+2\,B\,a\,b\,e^2\,f^3-2\,A\,a\,b\,e\,f^4+6\,B\,a\,c\,d\,e\,f^3-4\,A\,a\,c\,d\,f^4-2\,B\,a\,c\,e^3\,f^2+2\,A\,a\,c\,e^2\,f^3+3\,B\,b^2\,d\,e\,f^3-2\,A\,b^2\,d\,f^4-B\,b^2\,e^3\,f^2+A\,b^2\,e^2\,f^3+4\,B\,b\,c\,d^2\,f^3-8\,B\,b\,c\,d\,e^2\,f^2+6\,A\,b\,c\,d\,e\,f^3+2\,B\,b\,c\,e^4\,f-2\,A\,b\,c\,e^3\,f^2-5\,B\,c^2\,d^2\,e\,f^2+2\,A\,c^2\,d^2\,f^3+5\,B\,c^2\,d\,e^3\,f-4\,A\,c^2\,d\,e^2\,f^2-B\,c^2\,e^5+A\,c^2\,e^4\,f\right )}{f^5\,\sqrt {4\,d\,f-e^2}} \]

input
int(((A + B*x)*(a + b*x + c*x^2)^2)/(d + e*x + f*x^2),x)
 
output
x^3*((A*c^2 + 2*B*b*c)/(3*f) - (B*c^2*e)/(3*f^2)) + x*((A*b^2 + 2*A*a*c + 
2*B*a*b)/f - (d*((A*c^2 + 2*B*b*c)/f - (B*c^2*e)/f^2))/f + (e*((e*((A*c^2 
+ 2*B*b*c)/f - (B*c^2*e)/f^2))/f - (B*b^2 + 2*A*b*c + 2*B*a*c)/f + (B*c^2* 
d)/f^2))/f) - x^2*((e*((A*c^2 + 2*B*b*c)/f - (B*c^2*e)/f^2))/(2*f) - (B*b^ 
2 + 2*A*b*c + 2*B*a*c)/(2*f) + (B*c^2*d)/(2*f^2)) - (log(d + e*x + f*x^2)* 
(B*c^2*e^6 - 4*B*a^2*d*f^5 - A*c^2*e^5*f - A*b^2*e^3*f^3 + B*a^2*e^2*f^4 + 
 4*B*b^2*d^2*f^4 + B*b^2*e^4*f^2 - 4*B*c^2*d^3*f^3 + 6*A*c^2*d*e^3*f^2 - 8 
*A*c^2*d^2*e*f^3 - 5*B*b^2*d*e^2*f^3 - 8*A*a*b*d*f^5 - 2*B*b*c*e^5*f + 13* 
B*c^2*d^2*e^2*f^2 + 2*A*a*b*e^2*f^4 - 2*A*a*c*e^3*f^3 + 8*A*b*c*d^2*f^4 - 
2*B*a*b*e^3*f^3 + 8*B*a*c*d^2*f^4 + 2*A*b*c*e^4*f^2 + 2*B*a*c*e^4*f^2 + 4* 
A*b^2*d*e*f^4 - 7*B*c^2*d*e^4*f - 10*A*b*c*d*e^2*f^3 - 10*B*a*c*d*e^2*f^3 
+ 12*B*b*c*d*e^3*f^2 - 16*B*b*c*d^2*e*f^3 + 8*A*a*c*d*e*f^4 + 8*B*a*b*d*e* 
f^4))/(2*(4*d*f^6 - e^2*f^5)) + (B*c^2*x^4)/(4*f) + (atan(e/(4*d*f - e^2)^ 
(1/2) + (2*f*x)/(4*d*f - e^2)^(1/2))*(2*A*a^2*f^5 - B*c^2*e^5 - 2*A*b^2*d* 
f^4 - B*a^2*e*f^4 + A*c^2*e^4*f + A*b^2*e^2*f^3 + 2*A*c^2*d^2*f^3 - B*b^2* 
e^3*f^2 - 4*A*c^2*d*e^2*f^2 - 5*B*c^2*d^2*e*f^2 - 2*A*a*b*e*f^4 - 4*A*a*c* 
d*f^4 - 4*B*a*b*d*f^4 + 2*B*b*c*e^4*f + 2*A*a*c*e^2*f^3 + 2*B*a*b*e^2*f^3 
- 2*A*b*c*e^3*f^2 - 2*B*a*c*e^3*f^2 + 4*B*b*c*d^2*f^3 + 3*B*b^2*d*e*f^3 + 
5*B*c^2*d*e^3*f - 8*B*b*c*d*e^2*f^2 + 6*A*b*c*d*e*f^3 + 6*B*a*c*d*e*f^3))/ 
(f^5*(4*d*f - e^2)^(1/2))